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ABSTRACT 

!  In  many  applications  one  meets  systems  of  differential  equations  which 
consist  of  first-order  hyperbolic  and  seoond-order  parabolic  subsystems  which 
are  nonlinear ly  coupled.  These  arise,  for  instance,  in  the  modeling  of  motion 
of  a  compressible,  viscous  heat  conducting  fluid,  in  radiation  hydrodynamics, 
and  in  the  theory  of  motion  of  viscoelastic  materials .  The  relevant  equations 
are  presented  he*e±rrr:J_ 

The  results  of  this  work  are  local  time  existence  and  uniqueness  theorems 
for  initial-boundary  value  problems,  including  cases  with  free  boundaries,  for 
such  systems.  The  results  given  are  for  the  case  of  one  space  dimension.  The 
methods  used  involve  introducing  appropriate  variables,  the  method  of 
iteration,  a  priori  estimation  and  fixed  point  theorems.  -(  ■ 
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BOUNDARY  VALUE  PROBLEMS  AND  FREE  BOUNDARY  PROBLEMS 
FOR  QUASI LI NEAR  HYPERBOLIC -PARABOLIC  COUPLED  SYSTEMS 

Ta-tsien  Li*,  Wen-tzu  Yu**,  and  We-shi  Shen*** 

i  •  ■  :.TrOLiucrroN. 

'.eiij  wc  handle  with  so-called  quasi  linear  hyperbolic-parabolic  coupled  systems  which 
occur  in  applications.  Roughly  speaking,  in  this  kind  of  system  a  part  of  equations 
formulate  a  first  order  quasil inear  hyperbolic  system  with  respect  to  certain  unknown 
functions  x  «  {uj,***,u  ),  another  part  of  equations  a  second  order  quasilinear  parabolic 

system  with  respect  to  the  remainder  of  unknown  functions  v  *  (v4,...,v  },  and  these  two 

I  m 

.'arts  arc  non linearly  coupled  each  other.  For  instance,  the  system  of  motion  for  a 
compressible  viscous,  heat-conductive  fluid^1^,  the  system  of  radiation  hydrodynamics^^, 
r.he  system  of  motion  of  viscoelastic  materials ^ etc.  are  of  this  kind* 

The  initial  value  problem  with  smooth  initial  data  has  been  studied  by  several 
r.  j‘ ,hor3.  For  example,  for  the  system  of  motion  for  a  compressible  viscous,  heat-conductive 
fluid  ii.  1-dimensional  case,  J.  Nash  ^  and  N.  Itaya^^  have  proved  the  existence  and 
V'-v  uniqueness  c>f  the  local  smooth  solution.  Recently*  A«  Katsumura  and  T*  Ntshida^ 

;  j  wen  proved  the  corresponding  global  existence  theorem  for  the  small  initial  data. 

Fm  the  quasilinear  hyperbol ic-parabol tc  coupled  system,  the  boundary  value  problems, 
OfSpt-t illy  the  free  boundary  problems  are  more  important  in  applications,  because  the 

is  conce  rned  with  determining  the  corresponding  discontinuous  solution  which  can 
•  j'T.oribo,  *.  instance,  the  behaviour  of  a  fluid  containing  a  radiation  shock  in  radiation 
hydrodynamics.  But  for  the  boundary  value  problems,  especially  for  the  free  boundary 
.;'n4;  ‘-.as,  can  only  find  certain  results  in  some  special  cases  even  for  one-dimensional 
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case  (for  instance,  A.  Tani 


has  discussed  the  mixed  initial-boundary  value  problems  for 


(8) 

the  system  of  compressible  viscous,  heat-conductive  fluids  in  a  cy Under ical  domain  with  a 

♦ 

special  Dirichlet  boundary  condition:  the  velocity  u  ■  0  and  the  absolute  temperature 
T  ■  Tj(t,x)j  A.  V.  Kashikhov  and  V,  V.  Shelukhin^  have  considered  the  corresponding 
one-dimensional  initial-boundary  value  j  blem  with  the  boundary  data: 


u(t,0)  «=  u(t,1)  -  Tx(t,0)  *  Tx(t,1)  -0,  t  >  0> 

Moreover,  A  Tani^15^  has  also  studied  a  free  boundary  value  problem  for  compressible 
viscous  fluid  motion  etc.).  So  it  is  worthwhile  to  carry  out  a  systematic  research  on  this 
subject. 

In  what  f ol lows  we  shall  concentrate  our  attention  on  the  boundary  value  problems  and 
the  free  boundary  problems  for  the  following  general  types  of  quasilinear  hyperbolic- 
parabolic  coupled  systems  in  one-dimensional  case: 


n  3u 

I  CA  (t.x.u.y)^  ♦  \t 
J-1  J 


(t 


,x,u,v,vx 


Pt(t,x,u,v,vx)  (i 


dv 

Tt 


a(t,x,u,ux,v,v  ) 


b(t,x,u,ux,v,vx>  , 


1,. ..,n)  (1.1) 


(1.2) 


T 

where  v  -  (v v  )  is  a  vector  function  and  a  is  a  diagonal  matrix: 
i  m 

a  *  diag( a  , . . . ,a  )•  On  the  domain  under  consideration,  we  suppose  that 
I  m 

det|Ctj|  *  0 

and 

a^  )  0  ( 1  “  1 , . . • ,m)  . 

T 

In  this  system,  (1.1)  is  hyperbolic  with  respect  to  u  »  (u^,..*,un)  (under  the 

characteristic  form  with  the  characteristic  directions  ■  A,  (4  *  1,...,n)),  (1.2)  is 

at  * 

parabolic  with  respect  to  v  and  (1.1),  (1.2)  are  non linearly  coupled  each  other. 
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Type  XI. 


,u,v,vx>  **>- 

;1(t,x,u<v)(|^  +  Xj(t,x,u,v,v^)  J-j)  +  Mt(t,x,u,w,»x)  (»  - 


(1.3) 


3v 

Tt 


-  a(t,X,U,V,V  ) 


afv 

3x2 


b(t,x,u,v,vx) 


( 1.4) 


in  which  the  coefficients  a  and  b  don't  depend  on  ux,  but  on  the  right-hand  side  of 

(1.3)  there  is  an  additional  term  which  denotes  the  directional  derivative  of  v  along  the 

characteristic  direction  “  X,. 

dt  1 

Our  goal  is  to  discuss  various  kinds  of  boundary  value  problems  and  of  free  boundary 
problems  for  these  systems  in  a  class  of  smooth  functions  or  piecewise  smooth  functions  and 
give  a  condition  of  local  solvability  in  order  to  obtain  the  corresponding  existence  and 
uniqueness  theorem.  The  results  obtained  by  us  can  be  applied  to  many  practical  cases  and 
imply  an  affirmative  answer  for  a  conjecture  given  by  C.  M.  Dafermos^3'  about  the 
Incomplete  parabolic  damping. 


2.  EXAMPLES. 

1.  System  of  motion  for  a  compressible  viscous,  heat-conductive  fluid. 
In  one-dimenelonal  case  the  system  can  be  written  as  follows 


3p  . 

TZ  + 


3u 

re 


•f 


;  fc  ("  fe)  • ;  tr  * 1  ■ 


3t 

T[  * 


3T 


1 _  3  f.  3Ti  |i  f3u>2  '  P 

HZ  TrJ  ♦  7=5-  It:;  ♦  5“ 


PS 


P  3u 


(2.1) 

(2.2) 

(2.3) 


where  t  i  time,  x  s  spatial  coordinate,  P  :  density  (p  >  0),  u  i  velocity. 
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p  i  pressure,  T  <  absolute  temperature  (T  >  0) ,  St  entropy,  W  t  coefficient  of 
viscosity  (M  >  0),  X  t  coefficient  of  heat  conduction  (X  >  0),  ft  outer  force  which  is 
a  given  function  of  (t,x),  and  p,S,W,X  are  given  functions  of  (P,T). 

It  is  easy  to  see  that  in  this  coupled  system  (2.1)  is  a  single  first  order 
(hyperbolic)  equation  for  P,  (2.2),  (2.3)  is  a  second  order  parabolic  system  for  (u,T). 
So,  this  system  is  of  the  following  fora  of  quaailinear  hyperbolic-parabolic  coupled 
systems  t 


n  3u  j  3u 

l  etj(t,x,u,v)(^J-  +  Xt(t,x,u,v,-j^)  -  Ut(t,x,u,v,-5j)  (i  -  1,...,n)  , 


(*) 


-  a(t,x,u,-j^,v)  — j  -  b(t,x,u,‘j~,v,-jj)  , 
3x* 


in  which  Xt  and  U ^  are  affine  functions  of  r  ■  a  doesn't  depend  on  35  b 

does.  Obviously,  {*)  is  a  special  case  of  the  system  of  type  (I). 

2.  System  of  radiation  hydrodynaada. 

In  order  to  determine  the  motion  of  a  fluid  with  very  high  temperature  we  have  to 
consider  the  hydrodynamics  in  the  presence  of  a  radiation  field.  For  the  one-dimensional 
unsteady  flow,  under  the  diffusion  approximation  the  corresponding  system  of  radiation 
hydrodynamics  can  be  written  in  Lagrangian  representation  as  the  following  conservation 
law: 


3t  3u 

7?  '  5? 


0  , 


3u  i(p  *  pv* 

+  — 55 - 


0  , 


3(« 


2 

V 

Tt 


%  ♦  IF  )  3{u(p  +  pj 


DP 


3e 

Vi 

XT' 


"55 


in  which 

T  «  ~  :  specific  volume, 
p  -  RPT  :  pressure,  R  -  constant  >  0, 
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py  »  |  T3  :  light  pressure,  a  3  ,  c  s  light  speed,  0  *  constant  >  0, 

RT 

e  *  y— j*  :  inner  energy,  Y  :  adiabatic  exponent, 

4 

Ey  “  3pv  -  aT  :  radiation  energy, 
tC  .  0 

D  *  ,  *  ■  AT  (A,  a  >  0  constants):  Rosseland  mean  free  path. 

Taking  (u,p,T)  as  unknown  functions,  the  system  can  be  written  as 


'Sf  (g  ♦  *  P(fe  ♦  ■  -p(«p  *  ^  fT3)  g  . 

✓SF  C|H  -  p/5?  |E)  -  p(|a  -  p/i?  g]  -  P(RP  ♦  ii  §  T3)  g  . 


r  R  160  T3-,  3t  16AO  _ ,3+a  32T  16A0  _3+a  3p  3t 

^  +  —  r]  sF  •  ~  W  ^2--3~T  j:v; 


16AO 

3 


(3  +  a)pT2+“(g)2 


1  3F 


D  . 


(2.4) 


(2.5) 


(2.6) 


It  is  clear  that  (2.6)  la  a  single  second  order  parabolic  equation  for  T,  and  (2.4), 
(2.5)  is  a  first  order  quasllinear  hyperbolic  system  for  p  and  u  (with  the 
characteristic  directions  —  -  X,  .  ■  tP^OT) ,  then  (2. 4) -(2. 6)  is  also  a  special  case  of 

dt  1,4 

system  <*). 

3.  System  of  one-dimensional  viscoelastic  materials  of  the  rate  type 


+  p(u)  "  v 

t  x  x* 


(2.7) 


and  system  of  one-dimensional  thermovlscoelastlc  materials 


ut  -  vx  '  0  ' 


v  +  p(u,8)  ■  v 

t  X  XX 


(2.8) 


(e(u,9)  +  ♦  ip(u,  8)  v]  -  [nr  ]  -  8  (el  >  0) 

2  t  X  XX  XX  9 
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are  both  a  single  first  order  (hyperbolic)  equation  for  u  coupled  by  a  parabolic  equation 
or  system  respectively. 


System  of  one-dimensional  thermoelastic  materials 


0  , 


v  +  p(u,9)x  -  0  ,  (p^  <  0)  , 


[e(u,S)  ♦  |-]t  +  fp(u,8)v]x  -  8xj(  («£  >  0) 


(2.9) 


is  a  hyperbolic  system  for  u  and  v  coupled  by  a  single  parabolic  equation  for  0. 
These  systems  are  of  the  form  (*). 

System  of  a  model  of  nerve  impulse  propagation 


ufc  -  ♦  F0(u,w)  , 

wfc  -  G(u,w) 


(2.10) 


and  system  of  reaction-diffusion 


3u  3  u 


“  — r  +  f(u,v)  , 

3*2 


(2.11) 


3v  . 

-  g(u,v) 


are  obviously  of  the  form  {*). 

5.  Moreover,  certain  higher  order  aquations  can  be  also  reduced  to  a  hyperbolic- 
parabolic  coupled  system,  for  instance,  we  consider  the  following  problem  (see  J.  M. 
Greenberg^10'’  Greenl>er9»  R<  C.  MacCaray  and  V.  J.  Misel^1^,  also  see  J.  L. 

Lions1121  ): 

$-<&$-**&•••  •«■«•-  ••• 


u( 0,t)  -  u(  1  ,t)  -  0  , 

3u 

u ( x , 0 )  -  %(*>«  (**0)  ■  u,(x) 


(2.12) 
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Taking 


v 


3u 

it- 


3u 

w  ■ 


as  new  unknown  functions#  this  problem  is  equivalent  to  the  following  one 

3w 

t: 


3u  3w  3v  3v  .  3w  .  32v  . 

'  «i  tr"  t r>  ist  -  E(w)  ttt  “  *  — 7  “  0  , 

3x2 


Jt  *  'St  ”  Ti'  7t 

v(0,t)  •  v( 1,t)  -  0  , 


(2.13) 


u(x,0)  «  u0(*>»  w(x,0)  -  u^(x>,  v(x,0)  -  u1(x> 


in  which  the  first  two  equations  formulate  a  hyperbolic  system  for  u  and  w  and  the  last 
one  is  parabolic  for  v,  so  this  system  is  of  the  form  (*),  too. 

Now  we  shall  point  out  that  in  many  cases  by  means  of  adding  certain  new  unknown 
functions  some  problems  for  system  (*)  can  be  equivalently  reduced  to  a  corresponding 
problem  for  the  system  of  type  (II),  for  which  the  existence  and  uniqueness  theorem  seems 
somewhat  easier  to  prove. 

Example  1i  Consider  Cauchy  problem: 


n  3u  3u, 

I  Ctj(t,x,u,v)(^l  +  Xt(t,x,u,v,vx)  -  Mjl(t,x,u,v,vx),  (t 


,n)  ,  (2.14) 


3v 


32v 


x-  -  a(t,x,u,u  ,v)  — -  -  b(t,x,u,u  ,v,v  ) 
®t  *  Jx2  xx 


(2.15) 


t  •  0  :  u  »  *>(x) ,  v  -  ♦(x) 


(2.16) 


in  which  1^  and  *  1,...,n)  are  affine  functions  of  r 


w 


j 


(j  -  1,...,n)  , 


3v 

3J- 


Set 


differentiating  (2.14)  with  respect  to  x  and  using  equation  (2.15),  we  can  prove  that 
u,v  and  w  satisfy  the  following  Cauchy  problem 
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!Y7S5v 


V 


Ut(t,x,u,v,vx) 


3u . 


3u. 


Z1  C£j(t'*'u'v>(?p  ♦  At(t,x#u,v#Vx)  *5^i) 

n  3w  -■».  _  a  a 

I  +  \t(t,x,a,v,vx)  -  ?t(t,x,u,w,v)(^  +  \t 


3w, 


+  U.(t,x,u,w,v,v  ),  (t  -  1,...,n)  , 


(2.17) 


jv  j2y 

■J£  -  a(t,x,u,w,v)  — -  -  b(  t ,x,u,w, v,v^)  , 


t  »  0  1  u  «  *(x),  w  «  ^  '  (  X)  ,  V  «  9C) 


in  which 


Ct(t,x,u,w,v) 


n 

I  tt:)  w^ )/a(t,x,u,w,v) , 


(r  denotes  |^)  , 


U^t.x.u.w.v.r) 


3ut  n  3ut  3wt  jv 

sr +  Ji  ^  “k +  vr  x 


+ 


C,*,v 
*  t  x 


•4V 

are  determined  by  the  coefficients  and  (C  )  is  the  inverse  matrix  of  ( ) . 

Conversely,  if  (u,v,w)  is  the  solution  of  problem  (2.17),  then  we  can  prove  that 

9U 

(u,v)  i**  the  solution  of  the  original  problem  (2s  14) -(2. 16)  and  w  *  -r-* 

«x 

Example  2i  Taking 


8 


3w 

t: 


as  an  unknown  function,  problem  (2.13)  is  equivalent  to  the  following  problem  of  type  (II) 
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v?05»tfr 


3u 


3w  3v  3a 
v'  It  ’  Tx'  3t 


I  ^  -  E(w)s), 


3v 

sz 


E(w)  s 


v(0,t)  -  v( 1 , t )  -  0 


(2.18) 


u(x,0)  »  ufl(x),  w<x,0)  -  u’(x),  •(«»(>)  “  u;(x),  v(x,0)  -  u^  ( x)  . 

Here,  we  can  find  out  that  on  the  right-hand  side  of  the  third  equation,  there  is  a 

directional  derivative  of  v  along  the  characteristic  direction  •5^  -  0  and  that 

at 

/  3w  3s  \ 

a  •  X,  b  -  E(w)s  don’t  depend  on 

Hence,  in  order  to  explain  our  results  and  methods,  in  what  follows  we  shall  take  as 
an  example  the  second  initial-boundary  value  problem  for  the  system  of  type  (II).  all 
other  kinds  of  problems  (such  as  the  Cauchy  problem,  the  first  initial-boundary  value 
problem,  the  initial-boundary  value  problem  with  interface  etc.)  can  be  discussed  in  a 
similar  way  and  the  similar  results  for  the  system  of  type  (I)  hold  true,  too. 


3.  SECOND  IHITIAL-BOUNOARY  VALUE  PROBLEMS. 

On  a  rectangular  domain 

I 

R(«>  *  {( t,x) | 0  <t<«,  0<x«l)  I _ | _  (3.1) 

0  1  x 

we  consider  tl.s  second  initial-boundary  value  problem  for  the  system  of  type  (II): 


3v  . 32v 

JZ  -  a(t,*,u.v,vx)— 


b(t,x,u,v,vx)  . 


Without  loss  of  generality,  the  initial  conditions  may  be  written  as 

t-0:u*v«0. 


(3.3) 


(3.4) 
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Moreover  we  can  suppose  that 


a(0,x, 0,0,0)  s  1 

(Otherwise,  use  the  transformation  of  independent  variables 

dS 


/ 


and  that 


0  /a(0,5,0,0,0) 


b(0,x, 0,0,0)  =  0  , 


) 


(3.5) 


(3.6) 


Ctj(0,x,o,o)  =  S 


1, 

0, 


*  -  3  , 
*  *  j  ’ 


(3.7) 


(to  this  end,  it  is  sufficient  to  introduce  the  transformation  of  unknown  functions 

v  -  v  -  tb(0,x,0,0,0)  , 

_  n 

u,  -  I  C,,(0,x,0,0)u  )  . 
j-1  3  3 


Under  the  hypothesis  (3.7),  the  u^  [t  -  1,...,n)  are  called  the  diagonal  variables. 
The  boundary  conditio  ns  are  as  follows: 


on  x  »  1 , 


on  x  ■  0, 


u^  -  G^( t,u,v)  (r  -  1,...,hf  h  <  n)  , 

3v 

“  r+(t,u,v)  » 

u-  -  G-(t,u,v)  (s  =  k  +  1,...,n;  k  >  0)  , 
8  8 

3v 

3^  -  F_(t,u,v)  . 


(3.8) 


(3.9) 


(3.10) 


(3.11) 


Here  the  boundary  conditions  for  v  are  of  Neumann  type,  so  this  problem  Is  called  the 
second  initial-boundary  value  problem. 

We  assume  that  the  following  conditions  are  satisfied: 
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( 1 ) .  Conditions  of  orientablllty : 


A  (0,1, 0,0,0)  <  0,  A  (0,1,0,0,01  >  0 
r  s 


1 • fh 

h  ♦  1,. 


(3.12) 


A. (0,0, 0,0,0)  <  0, 


A. (0,0, 0,0,0)  >  0 


1, • • • ,te 

k  +  1 , . , . ,n 


(3.13) 


As  usual,  the  characteristic  directions  are  called  departing  characteristic  directions  on 
the  boundary,  if  as  long  as  time  increases,  they  go  towards  the  interior  of  the  domain. 


Thus,  on  the  boundary,  the  number  of  boundary  conditions  for  u  is  equal  to  the  number  of 
departing  characteristic  directions.  For  example,  on  x  •  1  the  number  of  boundary 
conditions  for  u  is  equal  to  h,  the  number  which  appears  in  (3.12). 

(2)  Conditions  of  compatibility! 

G  (0,0,0)  -  0,  G.(0,0,0)  -  0  (r  -  1,...,hf  <  -  k  +  1,...,n)  ,  (3.14) 

r  8 


3g  n  3g 

■sr -  (o,o,o)  +  l  <o,o,o)y  (o,i,o,o,o)  -  u  <0,1,0, o,o>  , 

3t  j.1  aUj  j  r 

(3.15) 

JG.  n  i& 

■j-2  (0,0,0)  +  [  T-5  <0,0, 0)u  (0,0, 0,0,0)  -  0.(0, 0,0,0)  , 

3t  j-1  3uj  3  * 

(r  »  1,...,ht  s  -  k  +  1,...,n)  , 


Ft<0,0,0)  -  0  .  (3.16) 

(3)  Conditions  of  smoothness!  the  coefficients  of  the  system  and  the  boundary  conditions 
are  suitably  smooth.  For  simplicity,  we  omit  the  detail  here. 
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By  Mans  of  certain  a  priori  estimations  for  the  solutions  of  the  heat  equation  and  of 
the  linear  hyperbolic  system,  using  an  iteration  method  and  the  Leray-Schauder  fixed  point 
theorem,  we  have  proved  the  following 

Theorem:  Under  the  preceding  hypotheses,  suppose  further  that  the  following 
conditions  are  satisfied: 


3g 

det|4rr<  "  *  0  <r'r*  "  1«***»h)  , 

3(3- 

det|4..,  -  (0,0, 0)|  *  0  <s,s*  -  It  +  1,...,n)  , 

ss  u*( 


(3.17) 


i.e.  the  boundary  conditions  may  be  rewritten  as 


on  x  ■  1 , 


ur  »  Hr(t,ua,v)  (r  -  1,...,h>  s  «  h  +■  1,...,n) 
3v 

5^  ■  V*.u.v)  » 


(3.18) 


then, 

R(4) 


on  x 


0, 


« I *  *  Hj(t,u«,v) 


(r  ■  1,...,k>  s  »  k  +•  1,...,n) 


(3.19) 


-  F_ (t ,u, v)  . 

the  second  initial-boundary  value  problem  admits  a  unique  local  classical  solution  on 
where  4  >  0  is  suitably  small. 


4.  IDEAS  OF  THE  PROOF, 

1.  A  priori  estimations  for  the  solutions  of  the  second  initial-boundary  value 
problem  of  heat  equations: 


Jv 


a  v 

— -  +  b(  t,x)  , 
3xZ 


t  -  0  :  v  -  0  ,  (4.1) 

*  “  0  1  *  ^1*t)»  x  -  1  i  • 

a 

1  2'° 
Suppose  that  on  the  domain  R(4^),  ^(t)  €  C  ,  “0  (i  -  1,2)  and  b(t,x)  e  C 


(0  <  a  <  1),  where 


Holder  space  of  functions  f  such  that  f  la  Holder  continuous  with  respect 


c6'“- 

to  t  and  to  x  with  the  exponents  8  and  a  respectively  (0  <  a,  0  <  1), 
then  It  Is  well  known^13'  that  problem  (4.1)  admits  a  unique  classical  solution  v  on 
R(4  )  with 

t  1  t 

V(t,x)  -  /  /  N(t,X;T,£)b<T,?)d5dT  +  /  N<t,X(T,1)*  (T)dT 

0  0  0 


-  /  N(t,x>T,0)*>  (T)dt  , 

0 


£  (t.x)  -  /  /’  b(T.«)d«dT  +  /"  a-N(tr.|’Tf1->  ^(T)dT 


0  0 


-  /*  ,i(T,dT  . 


(4.2) 


*  1  32n 


—7  (t.x:  -  /  /  — I  (t,X,T,E)(b(T,e)  -  b(T,x))dWT  +  /  N(t,X>T,n;  (T)dT 

3x2  0  0  3x2  0  2 


-  /  N(t.X,T.O)  ^(DdT 


(t,x>  -  /  /  (t,xjX,e)(b<T,o  -  b(t,e>)dtdT  +  /  N(t,x»o,e>b(t,C)de 


t  1 


3n 


0  0 


in  which 


t  t 

♦  /  N(t,x»T,1)J  <T)dT  -  /  N(t,X|T,0)J  (T)dT  , 

0  2  0 


^1(T)  -  ^  ^1(T)  (1  -  1.2)  , 


N(t,xiT,5)  -  I  [Gg (t ,x» T,2n  ♦  t)  ♦  GQ(t,x;r,2n  -  5)) 

n»-** 


(4.3) 


is  the  Neumann  function  for  the  second  initial-boundary  value  problen  of  the  heat  equation 
and 

(x  -  sr 

(4.4) 


G0(t,x|T,5) 


.  4,t-T>  (t»T, 


2/*(t  -  I) 

is  the  fundamental  solution  of  the  heat  equation. 
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(0  <  a  <  1) ,  where 


_  — 2+a 

Moreover,  on  R<°0)  v  8  C 


-2+a  r _  3f  3f  3  f  _  3f  „  _ 

C  -  tf|f,  3^,  j  continuoue,  ■g^  e  C 

3x 


1+a  i  a 

2'  3f  32f  1'“, 

'^'17  c  1 


(4.5) 


On  R(4),v4,  0<4<4g,  introduce  the  following  normal 

•f •  “  sup  |f(t,x) |  , 

( t ,x)  e  r<4) 


H®[f] 


|f(t 


eup 


1,x)-f(tJ,x)  |  B 


fc  <t  ,x),(t,,x)  (t.,— t  I 

’e  R (4>2  1  2 


,  H  (f) 


|f<t,x  )-f(t,x  ) I 

sup  - - -  ,  (4.6) 


(t,X  ,),(t,X,)  lx  -x  I 

’e  R(«)2  1  2 


hV]  -  H2[f]  +  H“[f]  (0  <  a,  8  <  1) 


Using  the  preceding  expressions  end  the  property  of  the  fundamental  solution,  after 
a  long  calculation  we  have  obtained  the  following  three  a  priori  eatimations  on 
R(4)  V  4,  0  <  4  <  4.1 


1* .  Ivl  ---  Ivl  +  ll^l  <  C,(«1/2  Ibl  +  l<l)  or  Ivl  <  C,(!^  Ibl  +  <■*!)  , 
def  «x  l  * 


(4.7) 


2*.  Ivl 


1  def 


Ivl  +  +  .ifi,  +  ,>2  .  *■ 


ST  *  *^T'  +  V  W 


<  C^flbl  +  42H*(b]  +  d^yl^)  (Iwl^  *  lyl  +  l>l)  » 


(4.8) 


—  2 

3*.  ,V,2  def  ,V,1  +  Ht2  ‘fe  +  H<S[lt1  +  H°[^2J  *  V,bl  +  H“lbl  +  '*V  • 

3x 


(4.9) 


in  which  (i  -  1,2,3)  signify  constants  depending  only  on  4^. 

2.  A  priori  estimations  for  the  solutions  of  the  following  initial-boundary  value 
problem  of  first  order  linear  hyperbolic  systems; 
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*  •ii  t"*1” 


"■T7* Jfv 


1  dv  .  3v 

a  C  function,  ^  4>  ia  only  continuous.  But,  noticing  that  it  is  the 

directional  derivative  of  v  along  the  characteriatic  curve  •  A  (t,x),  we  can 
integrate  by  parte  this  tern  Cg['£~  +  A^  ’|~)  when  we  integrate  the  system  along  the 
characteristic  curve,  then  we  can  prove  as  usual  that  problem  (4.10)  admits  a  unique 

1  1+B  1  fl  Q 

clessical  solution  u  on  R(40>  with  u0C  oruec  -  {f|f  e  C  ,  e  CP*  } 

under  the  different  hypotheses  of  smoothness  respectively. 

Introduce  the  following  classes  of  functions! 


F0  ‘  W*  ?t<  V  Ar(t,1)'  A  » ( t ,  0 )  ^  * 

3?.,  35,  35,  85, 

r  i-  *1  *1  _  _ *  _ *  .  1  l 

*1  “  >5tj'  3t  '  3x  '  U'  3t  '  3x  '  V  det|5tjl'  ' 


3Xi  3W 

r2  “  ri  u  h T'  V  ST’  rTtTTT'  TOo7} 


(t,j  ■  1,...,m  r  “  1,...,hi  a  “  m  ♦  1,...,n) 
and  the  following  norms  of  functions  on  R(4)  ( 0  <  6  < 


-  «u.  +  i£,  +  ,£l 


•  3u  3n 

-  'o*  ♦  + 


8  au  0  3u 

,u,i+6  - ♦  ».«$ +  ».<&  ' 


,ui  i+8  -  ,u,i +  hM?'  ♦  +  • 


H.8m  -  H*Cfl  ♦  H*[fl  (4.  17) 

and  the  constant  c  >  0  will  be  suitably  chosen  later  on.  By  means  of  the  integral 
relations  satisfied  by  u(t,x)  and  by  -jj  respectively,  after  a  long  calculation  we 

have  established  three  a  priori  estimations  on  R(4)  V  >  S  >  0  as  follows! 


1» 


(4.18) 


lul  <  <1  +  K1 5)  < +  (Hq  +  K,4)  Ivl  ♦  15,411)1 

in  which  constant  K,  depends  only  on  the  norm  IT^I  on  R(  6^ )  and 
H  -  2  sup  U.(t,x)  |  . 

0  t-1 . .  * 

(t,x)  e  R(«fl) 

2*.  »ul*  <  (1  +  dg’e  +  K2i8)«H  +  (*0  +  K2 <5 ) <  1  +  ly»1)  (4.19) 

provided  that  1^ , M ^  are  Holder  continuous  with  respect  to  t  with  the  exponent  8,  in 

which  d  ■  min  (-X  (0,1) ,X_(0,0) },  K  depends  only  on  the  norm  IT  I  on  R(  4  ) 

0  1<r<h  r  s  0  00 

k+1<s<n 

and  K2  depends  only  on  the  non  IT^I  and  H “ [ ]  on  8(4^). 

3*.  «ul*+B  <  (1  ♦  2  d"1*  ♦  dj2C  ♦  It/)..®!*]  +  (K2  +  Kj«)(1  ♦  l»l  ♦  lvl,+(()  (4.20) 


provided  that  all  the  functions  in  r  are  Holder  continuous  with  respect  to  t  and  x 
with  the  exponent  8,  where  Kj  depends  only  on  ir  I  end  on  R( Sg) . 

3.  Introduce  the  following  sets  of  functions  on  R( 5) I 


!.<«>  -  {(u.v)  luec1,  V  a  c1,  0  c°,  u(0,x)  -  v(o,x)  -  0}  , 

3x 

1^_a 

1,(4)  -  {(u,v)|u  e  C  2 ,  v  C  c"2+a,  u(0,x)  ■  v(0,x)  -  0, 

(4.21) 

3u  3 

<0,x)  -  U^IO.x, 0,0,0) ,  ( 0,x>  •  of  , 
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1 


where  A^,Bi  d  *  0,1,2)  are  positive  constants  to  be  chosen  later  with 

A0  ‘  A1  <  A2 '  R0  *  81  ‘  V 

For  any  (u,v)  8  X 1 < ® ) ,  according  to  the  preceding  points  we  can  define  an  iterative 
operator  (u,v)  -  T(u,v)  by  means  of  the  following  linear  problem 


n  „  «u  _  _  »u 

I  Ctl(t,x,u,v)(^l  +  Xt(t,x,u,v,|^)^J-) 


C((t,x,u,v)  (|^  +  lt(t,x,u,v,-jj)|j)  +  wt(t,x,u,v,|j)  (t  -  1,...,n)  , 


~ ?  -  b(t,x,u,v,|j)  +  [a(t,x,u,v,'|^)  -  1 J— r  . 
3x2  x  **  3x2 


t  -  0i  u  -  0,  v  -  0  , 


x  »  1»  j  C  <t,1,u(t,1),v(t,1)u.  •  G  (t,u(t,1),v(t,1)) 
j-1  3  1 


♦  Z  (C_,<t,1,u<t,1>,v(t,1)  -  S  >u.(t,1)  S  *  <t),  r  -  1 , . . . ,h  , 


(tf1)  «  F+(t,u(t,1),v(t,1))  2  *+<t)  , 


x  -  Oi  t  C- .<t,0,u<t,0) ,v(t,0) )u.  -  G.(t,u(t,0) ,v(t,0)) 
j-1  83  3  8 

+  5  (?gj(t,0,J(t,0),v(t,0)>  -  S^u^t.O)  2  »4(t). 


S  -  k+ 1 , • . • , n  , 


9  v  >  — 

•JJ  (t,0)  «  F_(t,u(t,0)  ,v(t,0)  2  *_(t) 


For  the  time  being,  we  suppose  thst 


n  3g  n  3G> 

I  IjT1  (0,0, 0)|  <  1,  [  |x-i  (0,0, 0)1  <  1 

j-1  Uj  j-1  #UJ 


.  -  V. 

rrr* 


Then,  using  the  preceding  a  priori  estimations,  we  can  choose  a  small  constant  e  >  0  and 


constants  such  that  the  operator  T  maps  £(4)  into  itself,  if 

4  >  0  is  suitably  small.  Because  £(4)  is  a  nonempty  convex,  closed,  compact  subset  of 
the  Banach  space  1,(4)  provided  the  norm 

•(u,v)l,  -  lul(  +  lv<1  ♦  I — jl 

3x * 

and  T  Is  a  continuous  mapping  from  £(4)  into  itself  in  this  space,  acc->r*  -  .  Leray- 

Schauder  fixed  point  theorem  this  operator  (u,v)  “  T(u,v)  has  a  fixe?  or.'.-'  j,v-  which 
is  the  solution  of  the  original  quasilinear  problem  on  R(4).  The  u,iitv>jei>i»r ;•  ti  he 
solution  can  be  proved  as  usual  by  means  of  the  corresponding  a  priori  'one. 

4.  In  order  to  finish  the  proof,  we  have  to  point  out  that  tie  ocatrae'  ton  condition 
(4.23)  can  be  realized  under  hypothesis  (3.17).  In  fact,  under  this  hypothesis  the 
boundary  conditions  can  be  written  as  (3. IB),  (3.19).  Then,  introducing  a  transformation 
of  unknown  functions 

ut  -  I1(x)u1  (i  -  1,...,n> 

and  multiplying  the  l-th  equation  of  (3.2)  by  I ^ ( x) ,  where 

Ii(x)  -  a^x  ♦  bjd  -  x) 

with  ag  “  bj  -  1,  af b^  >  0  small  enough,  it  is  easy  to  see  that  the  problem  for 
(u,v)  satisfies  the  corresponding  contraction  condition  (4.23)  as  well  as  all  the 
hypothesis  of  the  theorem,  so  the  theorem  is  proved. 


where  x  *  x(t)  is  an  unknown  boundary  curve,  we  consider  the  following  second  initial* 
boundary  value  problem  with  free  boundaries  (for  simplicity,  called  the  second  free 
boundary  problem) i 


3u, 


3u^ 


I  C^U.x.u.vK^1  +  Xt(t,x,u,v,vx)-j^-)  •  5t(t,x,u,v)  ♦  At(t,x,u,v,vx)|j) 


♦  U.(t,x,u,v,v  )  (t  »  , 


b(t,x,u,v,v  )  , 


and 


iv  .  3^v 

^  -  a(t,*<U,v.vx)— 

«x 


t  -  0 1  u  -  v  •  0  , 


x  -  0:  u-  •  G-(t,u,v)  (a  »  k  +  1,...,n>  k  >  0)  , 

3v 

5^-  “  F_(t,u,v)  , 


x  -  x(t):  u^  “  Gx(t,x,u,v)  (r  *  1,...,h>  h  <  n)  , 

3v 

-  F+(t,x,u,v) 


-  D(t,x,u,v,v  ),  x(0)  -  1 

at  X 


(5.1) 


which  la  an  ordinary  equation  to  determine  the  free  boundary  x  -  x(t)  in  the  procedure  of 
solution. 

This  kind  of  problem  can  be  met  in  the  motion  of  the  fluid  with  radiation  shocks. 

He  assume  once  more  that  the  corresponding  conditions  of  orlentability,  of 
compatibiity  and  of  smoothness  hold  true.  For  instance,  we  assume 


0, 1 ,0,0,0)  <  D(0, 1 ,0,0,0) ,  1^(0, 1,0, 0,0)  >  D(0,1,0,0,0>  , 

1.(0, 0,0, 0,0)  <  0,  X;(0, 0,0, 0,0)  >  0 


(5.2) 


(r  -  1,...,h,  s  »  h  ♦  1,.,«,ni  r  -  1,...,k,  s  -  k  +  1,...,n)  . 

Here  the  essential  difficulty  consists  In  the  presence  of  the  free  boundary  curve,  but. 


using  the  transformation  of  independent  variables 

t  -  t,  x 

the  domain  R(4)  is  reduced  to  the  domain 


x(t)  ' 
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(5.3) 


S(S>  -  (( t,x) 1 0  <t<9,  0<x*O 

n,.,  —tau..  «-  «.  th.  Mlllu 

”  x“'  “°h  ““  lh*'  “•  ■*•»««.  °<  <•.«.  ~t  u  „  „y.  „ 

.  UUMUa  »■!«•  W  1«  „  ,„1W.  .M, 

again  replaced  by  (t,x))i 


!,  5*j,t*x|u'v>^+  -  et(t,xiu.v>(|i+  xt,t#x|u,v>|v} 


+  Mjtt.xlu.v)  (*  ,  # 

Jv  j2 

^  -  a(t,x|u,v)  -  b(t,x|u,v)  . 

*x 

t«Oi  u  *  v  -  0  ,  (5.4) 

*  ’  1«  “r  *  Gr(t|u,»)  , 

3v 

^  -  Vt|u,v)  » 

X  -  O:  u.  -  G.(t|u,v)  , 

3v 

■5J  -  P_(t|u,v)  , 

where 


-21- 


I 


CJjj(t,x|u,v) 

Xt(t,x|u,v) 

?j(t,x|u,v) 

Pt(t,x|u,V) 

a( t,x|u,v) 

b(  t,x|u,v) 

G^( t | u,  v) 
Gg(t|u,v) 
F+(t|u,v) 
r  (t|u,v) 


«  Cjj ( t,x( t) x,u, v)  , 

»  (Xg(t,X(t)X,U,V,  *  x(t) )  -  x’(t)x)/x(t) 


ct(t,x(t)x,u,v)  , 


Mt(t,x(t)x,u,v,  £  •  . 


»(t,x(t)x,u,v,  -jj  •  jjrjjO/x  <*>  • 


.  r  1  \  .  x' (t)x  3v 

b(t,x{t)x,u,v,  ^  ♦  "^tT  K  ’ 


Gr(t,x(t) ,u,v)  , 


-  G-(t,u.v)  , 

«  r+(t,x(t) ,u,v) *x(t)  < 
•  P  (t,u,v) -x(t) 


(5.5) 


(5.6) 


and  x  -  x(t)  la  defined  by 

D(t,x(t),u(t,1),v(t,1),  jj  (5,7) 

x(0)  -  1  . 


For  the  aecond  Initial-boundary  value  problem  in  functional  fora  we  can  prove  that  the 
altuation  la  eirailar  to  the  aecond  initial-boundary  value  problem,  then  we  can  obtain  the 
corresponding  condition  of  solvability  for  the  original  second  free  boundary  problem  ae 
follows t 


3g. 

detl -  ■jjjS.  (0,0, 0)1  *  0,  (s,s  -  k  ♦  1,...,n) 


3G 

det|«  -  -  w— £  (0,1, 0,0)|  *  0  (r,r  -  1,...,h)  , 
rr  <>u- 


(5.8) 


i.e.  the  boundary  conditions  can  be  written  as 

x  -  0*  u-  -  H-(t,u~,v)  (r  -  8  -  k  ♦  1,...,n)  , 

a  s  r 

3v 

■jj  -  F_(t,u,v)  , 

x  »  It  ur  “  Hr<t#x«ua»v>  (r  "  a  -  h  +  1,...,n) 

3V 

■jj  -  F+(t,x,u,v)  . 


(5.9) 


6.  VARIOUS  REMARKS. 

1.  In  the  case  where  the  given  boundary  x  -  x(t)  is  the  k-th  characteristic  curve, 
we  can  consider  x  *■  x(t)  as  a  free  boundary  with  the  condition 

-  *k(t,x,u,v,  |j),  x(0>  -  1  .  (6.1) 

Using  the  preceding  tranaformation  we  obtain  again  a  second  nixed  initial -boundary  value 
problem  in  functional  form. 

2.  In  a  similar  way  we  can  also  solve  the  following  problems! 

I*,  the  Cauchy  problem) 

2*.  the  first  mixed  problem  with  the  boundary  conditions ■ 
x  -  Oi  u-  «  G-(t,u,v)  (s  -  k  +  1,...,n) 
v  «  F_(t,u) 

x  -  1i  u^  »  Gr<t,u,v),  (r  »  1,...,h)  (6.2) 

v  ■  F+(t,u) 

and  the  corresponding  first  free  boundary  problem. 

3* .  the  problem  with  the  interface  x  -  0  on  the  domain 

r<«)  -  R^fj)  u  R _(«> 

with 
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for  the  following  system 


RJi)  -  {( t ,x)  |  o  <  t  <  S,  -i  <  x  <  0} 
\(t)  -  {(t,x)  1 0  <  t  <  «,  0  <  x  <  1 } 


(6.3) 


-  4t)(t.x,u(±>.v(t))(^i+  l‘i)(t,x.u<i).v(±>,|^.}|^) 
u  *  *x  '*x  J 

*  «4±>(t.x,u(t),v(±,,|i-^)  (t  ,  f 


t  »  0:  «  v***  -  0  , 


x  -  0:  ul+)-  Gi(t,u(t\,(t))  s  -  k  +  1, 


on  ((*(«)  respectively. 


ur‘  -  Vt'u<*)'v<t>>  r  "  , 


'<+>  .  v(_)  +  f(t,u(t))  , 


3v 


(♦> 


a(t,u(t> ,v(t))  - *  glt.i*’  ,4**  > 


(-) 


(a  >  o>  » 


X  ■  il:  convenient  boundary  conditions 
with  the  following  hypotheses: 


(0,0, 0,0,0)  <  0,  W+><0,0, 0,0,0)  >0,  r  •  a  -  k  ♦  1 


^’’(O.O.O.O.O)  <  0,  ^-'(0,0,0, 0.0)  >  0,  r-1 . h,  .  -  h  ♦  1 . . 


(6.4) 


(6.5) 
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4*.  The  problem  with  free  interface  x  -  x(t)  can  be  similarly  discussed,  boo.  This 
*e  the  conditions  on  x  •  x(t)  are  the  fol lowing i 


x  -  x(t)  :  ul+)-  G-Jt.x.u^’.v1**), 


k  +  1,...,n 


u^_)  -  Gi(t,x,ut4,,v(±)),  r  -  1, 


v(+)  -  v<_)  +  T(t,x,n<4> )  , 


(6.6) 


1 


*v<+>  .  _  u(±)  u(±).9v(">  .  „  „<*>  (*>. 

y  “ —  ■  a(t,x,u  ,v  )  5  “  ♦  g(t,x,u  #v  ) 


(±> 


f  -  o(t.x.»'*V*>.!S-^).  x<0)-0 


(6.7) 


th  the  following  hypotheses: 

*j+>(0,0, 0,0,0)  <  0(0, 0,0, 0,0)  S  D(0 ) ,  lj+)(0,0,0,0,0)  >  D(0)  , 

X*_) (0.0,0, 0,0)  <  D(0 ) ,  (0,0, 0,0,0)  >  D(0)  (6.8) 

(r  »  1, ...,k;  s  -  k  +  1,...,n>  r  -  1,...,hi  s  -  h  +  1,...,n)  . 

:  the  systea  of  conservation  laws,  these  conditions  (6.6),  (6.7)  and  (6.8)  (with 
«  k  -  1)  can  he  obtained  from  the  corresponding  Ranklne-Hugonlot's  conditions  and  the 
rrespondlng  entropy  condition  respectively. 


Similarly,  we  can  consider  the  problem  with  the  characteristic  interface  x  -  x(t) 


7. 


APPLICATION  TO  A  CONJECTIVE  GIVEN  BY  C.  M.  DAFBRMOS 


The  conjective  given  by  C.  M.  Dafermos^  is  that  incomplete  parabolic  damping  can 
preserve  the  smoothness  of  smooth  initial  data  but  ia  incapable  of  smoothening  rough 
initial  data.  For  the  system  of  one-dimensional  viscoelastic  materials  of  the  rate  type, 
he  has  verified  that  this  conjecture  is  true. 

Now,  using  the  preceding  results  we  can  consider  this  conjecture  in  general  case, 
indeed,  the  system  with  incomplete  parabolic  damping  is  a  hyperbolic-parabolic  coupled 
system,  since  the  problem  with  the  free  interface  (or  characteristic  interface)  x  *  x(t) 
is  well-posed  (under  the  corresponding  conditions  of  solvability)#  if  the  (rough)  initial 
data  are  piecewise  smooth  with  a  discontinuity  at  the  origin  x  m  0,  satisfying  the 
corresponding  conditions  of  compatibility  (i.e,  corresponding  Rankine-Hugoniot's  conditions 
for  conservation  laws)  and  the  corresponding  conditions  of  orientability  (i,a. 
corresponding  entropy  condition) #  then  the  local  solution  is  also  piecewise  smooth  with  a 
discontinuity  on  x  -x(t),  because  the  corresponding  conditions  of  solvability  can  be 
checked  in  many  concrete  cases.  Thus,  incomplete  parabolic  damping  is  incapable  of 
smoothening  rough  initial  data,  that  is  to  say,  the  second  part  of  this  conjecture  is 
true.  On  the  other  hand,  for  smooth  initial  data#  according  to  the  proceeding  results#  the 
solution  remains  smooth  locally  in  time,  so  the  first  part  of  this  conjecture  is  true  at 
least  in  a  local  sense.  As  to  the  corresponding  global  existence  theorem#  we  have  to 
discuss  the  concrete  system  and  the  problem  seems  yet  open. 
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